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Abstract. The spot model has been developed by Bazant and co-workers to describe quasistatic granular 
flows. It assumes that granular flow is caused by the opposing flow of so-called spots of excess free volume, 
with spots moving along the slip lines of Mohr-Coulomb plasticity. The model is two-dimensional and has 
been successfully applied to a number of different geometries. In this paper we investigate whether the spot 
model in its simplest form can describe the wide shear zones observed in experiments and simulations of a 
Couette cell with split bottom. We give a general argument that is independent of the particular description 
of the stresses, but which shows that the present formulation of the spot model in which diffusion and 
drift terms are postulated to balance on length scales of order of the spot diameter, i.e. of order 3-5 grain 
diameters, is difficult to reconcile with the observed wide shear zones. We also discuss the implications 
for the spot model of co-axiality of the stress and strain rate tensors found in these wide shear flows, and 
point to possible extensions of the model that might allow one to account for the existence of wide shear 
zones. 
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1 Introduction 

It is well known that a relatively general theory for flow 
of granular media is still beyond reach. In part, this is 
due to the richness of granular flow phenomena [UEKE] 
- ranging from avalanche type behavior down inclined 
planes or granular surfaces to hopper discharges, granular 
gases, sheared flows in Couette cells or glassy type rhe- 
ology. Other impediments to the development of a gen- 
eral framework to describe granular flow are the compe- 
tition between static regions and flowing regions, and the 
fact that flow zones are typically not much wider than 
a few particle diameters, so that continuum theories are 
questionable. Even though a general theory is still lack- 
ing, some more limited approaches aiming at describing 
the phenomena in a particular limit or in some special 
case have been relatively successful [5JIMI]- O ne particu- 
lar recent proposal that appears to be quite successful for 
describing hopper flow and similar quasi-static flow prob- 
lems is the so-called spot model introduced by Bazant and 
co-workers [8ll9|IT0 l [TT ) . The central idea of this approach, 
illustrated in Fig. [TJ is that the slow flow of dense random 
packings can be understood in terms of the drift and diffu- 
sion of relatively lower density regions of about 3-5 grains 
in diameter, called spots. While this idea is not unlike ear- 
lier approaches incorporating cooperative rearrangement 
effects, including the Soft Glassy Rheology phenomenol- 
ogy [T^irTS] and the Shear Transformation Zone Theory 
[14,15,16,17], an attractive feature of the spot model is 




Fig. 1. Caricature 8 9 of a spot, a region where the density 
of grains is less than in the surrounding. When the spot moves 
upwards, the individual grains move downwards. 



that it is more specific, permitting calculation of flow prop- 
erties within a number of different experimental geome- 
tries and hence comparison to experiments or simulations. 
A second appealing aspect of the spot model is that it aims 
to merge the behavior at the scale of a few grains to the 
more "engineering" continuum-type approaches relating 
slow plastic flows to the stress fields (e.g. Mohr-Coulomb 
theory) . 

In this paper we will explore to what extent the spot 
model can be applied to the wide shear zones found in the 
so-called split-bottom geometries, i.e., cylindrical Couette 
cells in which the bottom consists of two rings, the inner 
of which rotates with the inner cylinder, the outer with 
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the outer cylinder [T8l. 19 , 20 , 2T] . The shear zones found 
in such cells are exceptionally wide, much wider than the 
shear zones of order 5-10 particle diameters found in many 
other flow geometries [2]. In addition, the shear zone bends 
inwards towards the inner cylinder as a function of height 
in a way which appears to be determined largely by the 
balance of torques [22] . 

A feature of the experimental data that made us ini- 
tially optimistic that the spot model in its simplest form 
could be applied to these wide shear zones is an empirical 
finding, namely that the shear rate as a function of radius 
follows an error function profile with an effective width 
that depends on height [18] . As the spot model leads to 
a diffusion-type equation for the density of spots, funda- 
mental solutions of the diffusion equation like Gaussians 
or error functions can be expected to emerge quite nat- 
urally for the spot density and the associated flow field. 
This intuitive observation motivated us to investigate the 
applicability of the spot model to these wide shear zone 
flows. 

Contrary to our expectations, our main finding is that, 
without significant modifications or extensions, it is hard 
to reconcile the spot model with the main experimental 
observations on the wide shear zones. To guide our discus- 
sion with a concrete example, we will first briefly discuss in 
section 3 the results of a simple-minded layer approxima- 
tion, which permits straightforward extension of the main 
features of the spot model in two dimensions to the three 
dimensional split-bottom geometry. This crude approxi- 
mation always predicts quite narrow shear zones which 
essentially go straight up, rather than bend inwards as 
is found in experiments and according to the torque bal- 
ance argument [22]. As we shall see, these shortcomings 
are, however, not due to the inadequacy of the layer ap- 
proximation underlying the implementation, but instead 
are intimately related to the basic structure of the drift- 
diffusion equation for spots. Indeed, in sections [4] and [5] 
we explore in generality which features of the spot model 
may be responsible for the incompatibility of the model 
with the wide shear zones. 

There are essentially two different conceptual ingredi- 
ents of the spot model which can be investigated sepa- 
rately: the more fundamental postulate that the flow can 
be analyzed in terms of a diffusion-type equation for the 
spot density, and the additional postulate that the drift 
vector in this diffusion equation can be calculated approx- 
imately from a Mohr-Coulomb type theory that predict 
stresses in materials at incipient yield. As we shall dis- 
cuss, without additional modifications to the model, both 
features appear to be problematic for these wide shear 
flows. On the one hand, the drift-diffusion equation can 
be shown to be incompatible with the main experimental 
features of the wide shear zones. At the same time the co- 
axiality of the principal axes of the quasistatic granular 
stress and strain rate tensors, which according to theo- 
retical arguments [24] and recent simulations [25] holds 
quite well in these wide quasi-static shear zones, is vio- 
lated by Mohr-Coulomb stresses in the two dimensional 
Couette setup where the principal directions of the strain 



rate tensor are fixed by symmetry considerations. In the 
outlook we will briefly mention some of the possible mod- 
ifications and extensions that may be required to describe 
wide shear zones within a picture of drift and diffusion of 
spots. 



2 Essentials of the spot model 

As stated above, the essential idea of the spot model is 
that flow in granular matter is mediated by the opposing 
movement of spots, regions of excess free volume, on the 
order of 3-5 grain diameters in size. The excess free volume 
associated with a single spot is less than the volume of one 
grain. When a spot moves in one direction, a net particle 
flow is caused in the opposite direction. The movement of 
the spots is postulated to be the result of a combination 
of drift and diffusion; this is called the "stochastic flow 
rule" . Together with the incompressibility of the flow this 
leads to the following central equation 

relating the granular velocity field u to the dimcnsionless 
spot density p. Here L is the spot size, At the time it takes 
a spot to move a distance L, and d the normalized drift 
direction vector. The "bare velocity" u is smeared out 
by spatial convolution with a "spot influence function" to 
give the final flow field. Nevertheless, the essential features 
of the flow are captured by u and the drift vector d, on 
which we therefore focus. Note that the time scale At sets 
the velocity scale, so that only one parameter, the spot 
size L, governs the balance between the drift and diffusion 
terms. As we shall see, this simple feature is intimately 
tied to the fact that the spot model tends to give rise to 
narrow shear bands of width L « 3-5 grain diameters. 

The above expression can be thought of as the consti- 
tutive equation of the spot model: it postulates that gran- 
ular flow results from the combined effect of a systematic 
drive (the drift term) and a random diffusion. The equa- 
tion also illustrates that there are two sides to the spot 
model. On the one hand, Eq. ([1]) captures the essential 
idea that quasistatic granular flow can be captured com- 
pletely in terms of the drift and diffusion of spots of re- 
duced local density. This idea can already be compared 
with known flow profiles if we know the drift direction 
vector d from simple symmetry arguments. E.g., in a two- 
dimensional hopper discharge it is argued [U[9] that the 
drift direction vector d points straight up, as a result of 
gravity. In this case, the equation for the spot density 
reduces to a simple diffusion equation, with the vertical 
direction playing a role analogous to time. We will use 
similar symmetry arguments below in section [4] to investi- 
gate the compatibility of the spot model with wide shear 
zones. 

The second side of the model is evident in nontriv- 
ial geometries, in which it is necessary to calculate the 
drift direction vector d explicitly. Calculation of d requires 
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Fig. 2. (left) Vector field for the drift direction vector d in one 
quarter of a cylindrical Couette cell cross section (inner and 
outer radii Vi„ and r ou t, respectively) driven with azimuthal 
velocity vg t drive at the inner boundary. For a Couette geometry 
in the absence of gravity, the drift direction vector is identical 
in every cross section. Note that while the drift has a large in- 
wards radial component, this will be compensated by an equal 
and opposite diffusion, resulting in only azimuthal flow in the 
steady state, (right) The azimuthal flow velocity as a function 
of radius on a logarithmic scale, illustrating the almost purely 
exponential fall-off of the velocity, consistent with for example 
Ref. [2]. 

a theory of stresses in the material; it is here that the 
spot model in its simple form makes contact with tra- 
ditional continuum methods from engineering. In order 
to deal with general 2d or quasi-2d cases, Kamrin and 
Bazant p~0|JTT| have proposed that d can be determined 
from the Mohr-Coulomb plasticity theory (MCP theory). 
MCP theory posits that, immediately prior to plastic fail- 
ure, granular materials are at incipient yield everywhere, 
i.e. the material is just about to fail, collapse or exhibit 
plastic flow along at least two lines through each point of 
the material, called slip lines. The direction of these slip 
lines is determined by the stress fields and the value of the 
Coulomb friction coefficient (i. Kamrin and Bazant assume 
that the stress tensor obtained from MCP remains valid 
in dense, quasistatic flows. They assert that spots per- 
form a biased random walk along the slip lines, the bias 
being provided by the net force on a locally fluidized mate- 
rial region. This then determines the effective drift vector 
d from the stresses in the granular medium. Note that 
this formalism is essentially two-dimensional, a property 
that it inherits from the MCP theory. As an example of 
the type of result that can be achieved with the current, 
two-dimensional, spot-model we present both d and the 
azimuthal velocity in a two-dimensional Taylor- Couette 
disk in Fig. [21 both results can be found in Ref. [TP] . 

3 Example: crude layer approximation for the 
split-bottom geometry 

In order to give an idea of the difficulties of reproducing 
wide shear zones within the spot model, we briefly sketch 
the main result of a crude and minimal extension to three 
dimensions that we have developed to study the steady 
state shear profiles in the split-bottom Couette geometry 



Fig. 3. (a) A schematic drawing of the split-bottom Couette 
geometry or Leiden geometry. The central part of the bottom 
rotates; the walls and the outer part of the bottom are sta- 
tionary, (b) A cross-section of the flow. The center and width 
of the shear zone are indicated. Note that it is relatively wide 
and moves inwards and widens with increasing height. Figure 
taken from Ref. [18] 

of Fig. [3] [T8 l ll9 y 20 p 21|. The shear bands observed in this 
system can exceed 50 grain diameters, much wider than 
what one typically finds. The location of the center of the 
shear band in this system is relatively accurately deter- 
mined by a simple torque minimization argument applied 
to an infinitely thin shear zone [22] . One would like the 
three-dimensional spot model to reproduce this and at the 
same time smear out these discontinuities to the relatively 
wide shear zones in the velocity field that are found ex- 
perimentally to be well fitted by an error function profile. 

There arc a number of ways to expand the two-dimen- 
sional spot model to a three-dimensional one. One option 
is to replace the essentially 2d MCP theory with a 3d 
plasticity theory. Although possible, 3d plasticity theo- 
ries come at the expense of a cumbersome mathematical 
framework. Below we will identify difficulties with the spot 
model inherent to both the stochastic flow rule and MCP. 
Though some plasticity theories may avoid the difficulties 
we identify with MCP, none can correct the issues associ- 
ated with applying the stochastic flow rule to wide shear 
zones. Thus, for simplicity, we illustrate a generalization to 
3d that builds on the 2d spot model incorporating MCP. 
We stress once more that we choose this simple and crude 
implementation here only for illustrative purposes — the 
model can certainly be made more realistic but the argu- 
ments of sections 4 and 5 show that this will not change 
the basic tendency of the model to give straight narrow 
shear zones, as long as one sticks to Eq. (JTJ. 

For the Couette geometry with symmetry along the 
axial direction, Kamrin and Bazant [10| have already ap- 
plied MCP. The corresponding d- vectors are indicated in 
Fig. [2] We build on these results by using a weak gravity 
approximation in which we consider the granular material 
to be composed of thin stacked slices. In each of these slices 
we assume that the two-dimensional picture holds, i.e. we 
assume plane-strain boundary conditions. This means we 
can solve the two-dimensional spot model, based on MCP 
theory, in each slice to find the two-dimensional drift vec- 
tor that would be there if this slice was the entire system 
23J. To this two-dimensional drift vector we then add a 
small vertical component due to gravity, which is constant 
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both within each slice and between slices. The drift vector 
is thus 

d = d|| (r, 6, z) + d±z. (2) 

Each slice is divided into two rings: an inner one that tends 
to rotate with the inner cylinder, and an outer one that 
tends to remain stationary in the lab frame. The bound- 
ary between the two rings is pinned at the split bottom 
but free to move with increasing height. Its position is 
determined sclf-consistently as we iterate upwards slice- 
by-slice. 

Using this approach, one is able to calculate the az- 
imuthal velocity in a split-bottom Taylor-Couette setup 
as a function of both height and radius. As can be seen 
from Fig. [4] these results do not match the experimental 
observations of Fig. [3] the shear band does not move in- 
wards for larger heights and also shows no widening - 
its width remains of order L, i.e. of the order of a few 
particle diameters. In the weak gravity limit, the fact that 
the shear bands remain of order L can actually easily be 
understood from the result of Fig. [5^: we already argued 
that as long as the d-vector points mostly in the radial 
direction, the spot model constitutive equation ([1]) pre- 
dicts that the shear band width will be of order L in 
the radial direction. We have explored various parame- 
ter regimes and other approximations, including a large 
gravity limit, but they always lead to essentially the same 
result: we invariably find relatively narrow shear bands 
that shoot straight up. Since it was shown by linger et 
al. [22] that the position of the shear band as a function 
of height can be obtained from a torque minimization ar- 
gument, we consider it likely that it is possible to get the 
position of the shear band from a more accurate contin- 
uum stress calculation than we have done here. However, 
we focus here on the width of the shear band and now 
proceed to show that this is a generic feature of using 
spot model expression ([I]) also in three dimensions, inde- 
pendent of which theory of stresses, MCP or otherwise, is 
used. 



4 Stress independent analysis 

For simplicity, let us consider a linear split-bottom shear 
cell. As sketched in Fig. El this is an infinitely long rect- 
angular container filled with grains. The bottom of the 
container is split lengthwise, so that the two halves can 
be moved relative to each other. Simulations [25] have 
confirmed that in this rectangular cell one obtains wide 
shear zones just as in the cylindrical Couette cell, but the 
geometry is somewhat easier to analyze due to the higher 
symmetry. If we denote the long dimension of the setup by 
y, the vertical dimension by z and the cross-channel hor- 
izontal dimension by x, symmetry dictates the flow must 
be in the y-direction and that physical quantities like the 
velocity or spot density cannot have any y-dependence. 

In the central vertical y-z-plane along which the two 
halves of the setup are sheared, both d x and d y must be 
zero due to symmetry, and since d is normalized d z = 1 
there. Neither in the experiments nor in the simulations is 
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Fig. 4. Results obtained in a three dimensional weak gravity 
approximation for the flow in a split-bottom Couette geome- 
try. Plotted are the azimuthal velocities as a function of the 
rescaled radius for three different heights. The shear band re- 
mains narrow and centered above the split in the bottom of 
the cell. 

there any sign that there is a nonzero velocity u z at this 
center line; if we therefore impose u z = at this center 
plane, we obtain from Eq. (fljl: 



which trivially leads to the exponential height dependence 

p = Aei\ (4) 

with A an arbitrary constant. This equation illustrates 
clearly the tendency of the naive extension of the spot 
model to 3d to lead to exponential profiles of width of or- 
der L. For the split-bottom geometry, such an exponential 
dependence can obviously not be correct: in the experi- 
ments [18, 19 ,20 ,21] as well as the simulations [25], wide 
shear zones at filling heights of h max = 50 grain diameters 
are studied. Even with L = 5 Eq. ((4]) leads to an overall 
height variation of the spot density by an unrcalistically 
large factor of order e 20 rj 5 • 10 8 — to put this in perspec- 
tive, note that one spot is thought to contain a "free" vol- 
ume of about a fifth of that of a grain [9]. Moreover, since 
the shear velocity in the y direction away from the center 
line is proportional to p, as there are no gradients in the 
y direction, such an exponential dependence would imply 
an exponential variation of the u y velocity with height z, 
unless d y is small and has a counterbalancing exponential 
height variation. Similar arguments apply if we analyze 
the extension of the profiles in the lateral (cross-channel) 
direction: for the profiles to be wide in the cross-channel 
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Fig. 5. A number of shear free sheets in a linear shear cell, 
the shear free sheet basis is indicated with 61,2,3- Figure taken 
from Ref. [25]. 

direction, we can assume d x ~ — ax with a <g; 1 (this 
Ansatz would lead to Gaussian variation of p in the lat- 
eral direction) but the height dependence would remain 
exponential and close to that of (j4|). Hence we conclude 
that independent of the precise "flow rule" that relates d 
to the local stresses in the granular medium, the extension 
of the spot model, based on treating all components of d 
on an equal footing (as in Eq. [T]) , is incompatible with the 
existence of wide shear zones. The reason for this is that 
the balance of the drift and diffusion term in Eq. |T]) is gov- 
erned by the single spot diameter length scale L — given 
that the spot size is hypothesized to be of order 3 — 5 grain 
diameters this almost inescapably leads to narrow shear 
bands of only up to ten grain diameters wide. 

5 Co-axiality and Shear Free Sheets 

The above discussion is independent of which particular 
flow rule is adopted for the connection between the stresses 
and the drift vector d. It is nevertheless of interest to go 
back and discuss the relation between the principal stress 
and shear directions on the basis of what we know from 
recent theory and simulations 24.25,26 , and to explore 
the implications for theories on granular flow. 

In the original formulation of the spot model, the so- 
called co-axiality flow rule is explicitly rejected and re- 
placed by a flow rule that builds on MCP theoy. Co- 
axiality means that the principal axes of the stress tensor 
and the strain rate tensor are co-axial, i.e., aligned. How- 
ever, according to the the analysis of wide shear flows of 
Dcpken et al. [21] , there are various reasons to believe that 
co-axiality is actually a crucial feature of these wide shear 
zones; later simulations confirmed this picture surprisingly 
well (21] . 

In order to explore what this implies, we need to ex- 
plain briefly the Shear Free Sheet (SFS) basis introduced 
by Depken and co-workers. A SFS is a surface of con- 
stant velocity; for the case of the rectangular split-bottom 
geometry the SFS's are sketched in Fig. [5] These sheets 
naturally define the orthogonal basis spanned by the unit 
vectors e\ and 82 in the SFS (with £2 is taken in the di- 
rection of flow) and £3 perpendicular to the sheet. Since 
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Fig. 6. (solid curve) The value of ip, the angle between the 
radial axis and the axis of minor principal stress, as a function 
of radius in a planar Couette setup. According to the results 
of Ref. [24] this should always be close to jtv (dashed line) or 
f 7r in a wide shear zone. 

by construction the velocity is constant within each sheet, 
the strain rate within each SFS is also zero. An easy calcu- 
lation then shows us that two of the principal directions of 
strain are at angles of 7r/4 relative to £2 and while the 
third principal direction is the same as e\. In line with the 
theoretical arguments [23] , recent simulations [2"M2"6"] con- 
firm that co-axiality holds in these wide shear zones, and 
hence that the orientation of the axes of principal stress 
are fixed by the SFS's. 

To illustrate how different the Mohr-Coulomb picture 
is from co-axiality in wide shear zones, let us return briefly 
to the approximation in which we think of the cylindri- 
cal Couette cell as being built up from slices of the two- 
dimensional (disk-like) Couette setup. Since the flow is 
azimuthal, the SFS's in each slice are just concentric cir- 
cles. Since the principal directions of strain rate are at 7r/4 
angles to these sheets, and since the shear zones bend only 
slightly inwards towards the inner radius with increasing 
height, they make an angle close to ±7r/4 with the az- 
imuthal direction. Due to co-axiality we know that this 
must also be true for the principal stress directions. This 
means that the angle between the radial axis and the mi- 
nor principal stress axis is in reality always close to jir 
or jTT depending on which of the two is the major and 
which is the minor principal stress direction. In the two- 
dimensional Mohr-Coulomb theory, however, the principal 
stress axes point in very different directions throughout 
most of the layer, as Fig. [6] illustrates. In other words, any 
plasticity theory that presumes a violation of co-axiality, 
docs not appear to be a viable starting point for the de- 
scription of wide shear zones. Any theory based on drift 
and diffusion of spots will have to build in, or self-consistently 
lead to, co-axiality in the wide shear zones. 



6 Outlook 

In putting our results into perspective we would like to 
stress that the spot model was not formulated specifically 
to apply to wide shear zones. In spite of this, the error 
function shear profiles found experimentally and the sue- 
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cess of the model in capturing Gaussian hopper flow pro- 
files led us to become optimistic that the model could cap- 
ture the wide shear zones as well. Nevertheless, our ana- 
lysis shows that present formulations of the spot model 
cannot capture the physics of the wide shear zones ob- 
served in split-bottom experiments and simulations. The 
main reason is that the model in its current form is based 
on the interplay and balance of diffusion and drift. As 
Eq. (TT]) clearly shows, this balance is governed by the sin- 
gle length scale L of order 3-5 grain diameters. As a result, 
in its present form the structure of the model is such that 
it leads to localization of the shear bands on this same 
scale L. 

More generally, our analysis brings up the question 
whether static stress considerations can be used to deter- 
mine the properties of wide shear zones: In Mohr-Coulomb 
theory and in the present spot model, one attempts to cal- 
culate the flow fields from static stress fields in conjunction 
with an incipient yield postulate, but the co-axiahty of the 
stress and strain rate tensor suggests instead a picture in 
which the flow "self-organizes" to a co-axial state. 

To be more specific: we have based our discussion on 
the central result, Eq. ([1]), for the connection between the 
drift vector d and the velocity field. In their discussion 
of the applicability of the spot model to other flow ge- 
ometries, Kamrin and Bazant [10] have independently ad- 
vanced the idea that spot drift must align with special 
surfaces (related to "slip line admissability" ) [27]. If in- 
deed we postulate that the spot drift is orthogonal to the 
z-direction in the central yz plane of Fig. [5] then clearly 
the main problem — an unrealistic exponential variation 
of the spot density p with height — dissolves. If d is not 
along the z-direction in the central yz plane, then by sym- 
metry it must be zero. Kamrin and Bazant in fact already 
interpreted a region with a null drift vector to be an in- 
dication that the spot mechanism is weak (a similar situ- 
ation occurs in plane shear and inclined plane flow [10] ). 
So from that perspective too, one unfortunately has to 
conclude that wide shear zones include physics that go 
beyond the spot approach. Other possible routes to ex- 
tending the spot idea to accommodate wide shear zones 
might be to introduce spot generation and annihilation 
terms [35J (i.e. regions of high plastic strain may act like a 
source of spots) ; to allow for evolution in spot size, thereby 
implicitly introducing another length scale; or to use an 
approach as suggested in Ref. [TO] that reconciles Bagnold 
type rheology with the physics of the stochastic flow rule. 
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